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ABSTRACT 




The stresses near the root of an F8U-3 airplane wing loaded in 
torsion were found using the matrix-force method of analysis. The wing 
was represented by a simplified, idealized model consisting of axial- 
load carrying flange bars and shear carrying cover plates. The theoret- 
ical stresses in the idealized model were found by using the IBM 360/67 
computer to solve for the redundant loads. These stresses show favor- 
able comparison with the experimental stresses determined from strain 
gages mounted on the cover skin of the wing. This study shows that the 
matrix-force method is applicable to thick-skinned wings. 
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INTRODUCTION 



The purpose of this study is to predict the stresses in an F8U-3 
airplane wing loaded in torsion using the matrix-force method of analysis. 
The stresses determined analytically are compared with the stresses ob- 
tained experimentally in order to determine the accuracy and consistency 
of the results. 

The primary structure of the F8U-3 wing center section consists of 
a simple multicell torque box, as shown in Figure 1. Spanwise bending 
is carried by the relatively thick skins which are fabricated from 
7079-T6 aluminum alloy plate. ' The skins are machined to varying thick- 
nesses both spanwise and chordwise in order to obtain an optimum distri- 
bution of material. The skins are continuous from the wing fold rib 
inboard to the wing centerline and from the front beam to the rear beam. 

It is this area, between wing fold ribs, port and starboard, and Front 
and Rear Beams, fore and aft, which defines the wing center section, 
the area of interest in this study. The sweep of the wing is 42® at 
the quarter chord. The wing thickness varies from 4% at the tip to 
5.01% at the root. It is this combination of a high degree of sweep 
and low thickness ratio which raises problems with structural analysis 
of the wing. Due to bending- tors ion interaction, or sweep-coupling, 
structural analysis by engineering beam theory methods runs into diffi- 
culty when analyzing stresses near the wing root. 

This thesis describes the matrix-force method of analysis for 
determining the stress distribution in such a wing. Among the advan- 
tages of this method are: 

(1) It results in a consistent set of internal stresses and external 
deflections . 
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FIG. 1 

F8U-3 WING CENTER SECTION 



12 



(2) The structure can be altered and re-analyzed with only minor 



changes in the computer data cards. 

(3) The ’^Bookkeeping” can be simplified and organized in a logical 
manner . 

The matrix-force method of analysis uses forces as the unknowns and 
treats the structure as an assembly of elastic components. The method, 
as described in Reference 1, is divided into three parts: 

A. Idealize the structure 

1. Replace the actual structure with a mathematical model. 

2. Make the model statically determinate by cutting redundant 
members . 

3. Write the equilibrium equations that determine the load distribu- 
tion in the structure 

B. Calculate and tabulate the flexibilities of the individual members 
comprising the structure. 

C. Perform matrix operations. 

In this report, the wing is ”cut” outboard of the intermediate 

rib (Y = 81.98), and it is assumed that the stress distribution here 
w 

can be accurately calculated using the engineering theory of torsion. 



assumed to be supported on four rigid supports located at A, B, C, 
and D in Figure 2. The aft supports are pins which can support fore 
and aft and vertical loads. The forward supports can support vertical 
loads only. The only loading condition considered in this study is 
S 3 nnmetrical wing torsion. From symmetry, the reactions at all four 
support points are determined by equilibrium as functions of the applied 
loads . 




dihedral of the wing is neglected, and the wing is 
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FIG. 2 

F8U-3 WING CENTER SECTION 
PLANFORM GEOMETRY a MEMBER LOCATION 



This study was conducted at the Naval Postraduate School, Monterey, 
California, during 1967. The author is indebted to Professor C. H. Kahr 
who conceived the need for the investigation and who rendered valuable 
assistance and encouragement throughout the project. The author expresses 
gratitude to Aeronautical Laboratory Technicians R. A. Besel and T. B. 
Dunton whose assistance and cooperation were essential to the completion 
of the study. 
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The model selected to represent the wing structure consisted of 
flat plates which carry only shear and cylindrical bars which carry 
only axial load. This model was selected for compatibility with the 
previous work done in Reference 2, Analysis of the model can be as 
accurate as the matrix-force method and arithmetic permit; however, 
the correlation between theoretical and experimental stresses will 
depend upon how closely the model represents the actual structure. 

This study is not an attempt to define the best model for the analysis, 
but to calculate stresses in what was tacitly assiimed to be an adequate 
model. Most of the examples of similar studies currently available in 
the literature deal with thin-skinned shells. There is little evidence 
on whether or not such models are suitable for the thick-skinned wing 
plates used in this study. An important factor in constructing the 
idealized model is the change of the wing geometry to model geometry. 

The model consists of bars and flat shear carrying plates. It is ex- 
tremely important to insure that these plates are, in fact, mathemat- 
ically flat. For example, the depths of the ribs at WS 25.095 and WS 
81.98 were such that the four points AFGH in Figure 1 did not lie in a 
plane. Therefore, panel AFGH is warped and has panel edge shears which 
have components perpendicular to plane AFGH. These shears are not 
accounted for in the equilibrium equations written for a flat panel 
and will produce errors in the solution. Reference 3 shows that for 
a similar wing such discrepancies yielded determinate member loads 
which were out of equilibrium as much as 5 to 10%. After altering 
the beam depths to make the panels flat, this trouble disappeared. The 
alteration required a change of less than 1% of each dimension. 
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After determining model geometry from the adjusted dimensions of 

the wing, the cross sectional area of the bars and the thickness of the 

shear panels were calculated. Each bar area of the model consists of 

the area of the beam flange plus one half the skin area between adjacent 

beams. The beam web areas were not included, using the methods of 

Reference 4. The effect of combining the axial load carrying skin with 

the beam flange to form the model's cylindrical bar is to remove Poisson's 

effect from the analysis. This means that any chordwise contraction of 

the cover skins is neglected. The wing skin thickness varies in the 

chordwise direction in section ABCD and in both spanwise and chordwise 

directions in section ABEF. An average skin thickness between beams 

was determined and used as the shear plate thickness. Using constant 

thickness shear plates in place of tapered skin simplifies the analysis 

but introduces possible sources of error. Model geometry and areas are 

listed in Tables I and II. Diagrams of cross sectional geometry are shown 

in Figure 3, 4, and 5. An exploded view of the outer section of the 

idealized model from Y = 25.095 to Y =81.98 is shown in Figure 6. 

w w 

An exploded view of the intermediate rib at Y^ = 81.98 is shown in 
Figure 7. 

2. Make the model statically determinate by cutting redundant members. 

The most complex beam that is truly determinate, so that the flange 
and shear loads can be determined from equilibrium equations, is the 
3-flange, 3 web beam shown in Figure 8. 




DETERMINATE BEAM 
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FIG. 3 

F8U-3 WING CENTER SECTION 
GEOMETRY OF CUT AT W.S. 25.095 
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F8U-3 WING CENTER SECTION 
GEOMETRY OF CUT AT W.S. 81.98 
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FIG. 5 

F8U-3 WING CENTER SECTION 
GEOMETRY OF CUT AT RIB (Y=0) 
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Fig. 7 

Exploded View of 
End Rib 



In this beam, loads can be obtained from equilibrium only, without con- 
sidering flange areas, material stiffness, or deformation of the structure 
In order to make the idealized model of this study a determinate structure 
it was necessary to cut the 5 interior webs and 11 of the flanges, re- 
sulting in the structure shown in Figure 9. 




■F/on^e. 27 
F/on^e 26 



Fig. 9 



CROSS SECTION OF DETERMINATE SYSTEM 



The number of redundancies arising from the geometry of a box structure 
has been given by Reference 4 and others as B + N - 4 where B is the 
number of longitudinal effective flanges continuous across each junction 
and N is the number of closed cells stiffened by ribs at tip end of bay. 
For the outer section of this model, B = 14 and N = b^which results in 
16 degrees of redundancy. This agrees with the 16 members which were 
cut to provide the determinate structure of Figure 7. 

3. Write the equilibrium equations that determine the load distri- 
bution in the structure. 

Considering the outer wing section from station 25.095 to 81.98, 
there are the following 49 unknown member loads. 

12 skin panels 

14 beam flanges 

7 beam webs 

10 rib flange loads at station 81.98 

6 rib webs at station 81.98 

49 total unknowns 
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The redundant forces come from the following members. 

1 F.B. , upper flange 

2 FR.I.B., both flanges 

2 F.I.B. , both flanges 

2 C.I.B. , both flanges 

2 A.I.B. , both flanges 

2 R.I.B. , both flanges 

5 webs 8-12 

16 redundant forces 

Thus, there are 49 - 16 or 33 equations of equilibrium required to 
characterize the model. These equations are: 

a) Six equations of equilibrium of free body of bay between stations 



Y = 25.095 and 81.98 



w 




(1) 


P = 0 

X 


(2) 


P = 0 

y 


(3) 


P = 0 
z 


(4) 


M = 0 

XX 


(5) 


M = 0 

yy 


(6) 


M = 0 
zz 



b) Eleven equations of equilibrium of the redundant beam flanges, elements 

16, 18, 20, 22, 24, 29, 31, 33, 35, 37, and 39. 

c) Ten equations of forces and moments, F = 0 , M =0 of the 

X ’ yy 

intermediate beam segments of the rib at Y =81.98 

w 

d) Six equations of equilibrium of the rib posts on the beam ends at 

Y^ = 81.98, elements 53 through 58. 

The addition of the center section Y = 0 to 25.095 can be thought 

w 

of as an extention of the horizontal flanges. The forces in the X and Z 



directions are reacted by the horizontal and vertical supports at Y^ = 
25.095, but the Y forces carry through to the centerline, Y^ = 0. 
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Although the center section is 16° redundant (the same as the out- 
board section), these are not new redundancies but merely components 
of the forces in the outer section. Further investigation shows that for 
symmetrical wing loading, the bending meoment, M , is constant from 

Y = 25.095 to Y = 25.095 and that torsional and vertical shear are 
w w 

zero. The axial load on each flange at Y = 25.095 is the same as the 

^ w 

load on each flange at Y = 25.095. Hence there is no shear on any 

w 

panel or web in the center section. (It was later shown that due to 
a slightly uns 3 nnmetrical load, there are, in fact, shearing stresses 
in the center section panels.) Since the addition of the center section 
results in no new redundancies and no new external loads, there can be 
no more than the 33 independent equilibrium equations already proposed 
for the outer section. Since there are more than 33 structural elements 
involved, there are several possible combinations of 33 independent 
equations. Those chosen were considered to be those most likely to 
have the larger stress concentrations. 

B. Calculate and tabulate the flexibilities of the individual members 
comprising the structure. 

As described in Section C, the internal strain energy of the wing 
may be written as: 



^£ir = Lo-, Qz ••• j w. 







( 1 ) 



The elements in the F matrix define the elasticity of the structure. 
If skin thicknesses or flange areas are changed, only the elements of 
the F matrix will be affected. If the geometry of the model is changed, 
the equations of equilibrium will also be affected. Multiplying out 
equation 1, 
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( 2 ) 



2 t. X 

£EU = t2a,a^ F,x 



i- 2a, a j i- •• ■ 



Each term of Equation 2 can be regarded as a small increment of strain 

2 

energy ^ t/ associated with the load a. or a, . a.. Since the load a. 

1 1 j 1 

may be applied to more than one structural element, the corresponding 

F.. term may be made up of contributions from several structural elements. 
11 

The F mattix must be a square symmetric matrix where F. . = F... 

Ji 

1. Flange Bars. 

The calculation of the F matrix elements for the flanges was done 
following the example of Reference 4. The energy in a uniform bar under 
a linearly varying axial force is 




2 






L. 



X 



dx- 



(3) 



where and <^^are end forces as shown in Figure 10. 
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Fig. 10 
FIANCE BAR 

Rewriting Equation 2 for this single element 

2 ^ ^ ^ 



2EU = t 2 % 2 



therefore 



^ a 



JL 

L 



X 

L 



dx 



f = — = f 

3A 



(4) 

(5) 

( 6 ) 
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and 



dx 






2/ 



'Sv 





( 7 ) 



In the case of linearly tapered members, the flexibilities may be 
obtained as functions of the end area ratios. Reference 1 shows the 
flexibility influence coefficients modified by a function, Q, of the 
end areas which was obtained from a graph of Q versus A. ^ 
area ratio. 



f = — (f> f - — ^ S 

For most of the flange elements, the generalized end loads, 
and were linear functions of the unkown member loads, or: 









( 8 ) 



(9) 



^ ^ ^ ■ ■ ' ( 5 *^ ( 10 ) 

where = the coefficient of the i*th member load applied to one end 
of the bar, 

and = the coefficient of the i'th member load applied to the other 
end of the bar. 

Now taking each force acting on the ends of the flange individually, 
that is : 

f (11) 

( 12 ) 

and substituting into Equation 4, we find that 

££(/ = -h p" ) (13) 
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( 14 ) 



and from Equation 1, for a single force, a^^ 

a£U = F„ 

Comparing Equations 13 and 14, 

Ffi — ^ 1 ^/ 

Similarly, for the other forces acting on the same flange, 

* -2°^* (3, ’Cl 

and for the cross product term, 

^ -f- ( 17 ) 

The forces on element 14 are shown in Figure 11. 

For example, the flexibility influence coefficients for element 14 were 
determined as follows: 






ft 

1 

_1 
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ll 




A^^/.3 


^ / y. j • 



Fig. 11 

FORCE DIAGRAM OF FLANGE BAR 

h ^ 'P,, = -iw - 



3^, 
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%, = -3.z 3/ - /o. 3 S3 

F =-Z0. ?33 (-/o. 9S3) fZx /0./67C-/0. 5S3)('o) +ZO. 533(^) “ Z 4 Z 9 > 

2 

^^2^"^o. 3Z3(~3.Z30 -h 20.323CoXs€e)P /0./<i7(-'0.9S3y:S66')=6S6.33S 
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For element 14, there are 4 more coefficients to be found; F 



^^ 21 , 21 ’ ^ 2 , 21 ’ ^^ 20 , 21 ' ^^ 2 , 2 



20 , 20 ’ 

coefficient for element 14 is 



added to other F^ « coefficients for elements 7, 39, 101, and 126 to 
make up the total F 2 2 flexibility influence coefficent found in the F 
matrix, 

2, Symmetrical Shear Panels 

Reference 3 gives the energy of the synimetrical trapezoidal 
panel in terms of the shear force on the shorter parallel side, as: 



2E\ 

so that 





/tv, X V 


} ’ 


" 1 






F 


i'u 
= F 




^ , n 


dci ^ 





may be written for F^ ^ of element 1 



F = 



\^4€. 



33S (/Z. 4c6) 



S/£ C. Zoo) 



■i(ooe7e) 







a - .8Z7a. 

C f 



(18) 



(19) 



Other F- ^ coefficients appear from elements 40 and 51 and add to pro- 
1 > ^ 

duce the final F- - found in the F Matrix. 
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The shear panel edge loads were determined by assigning a unit 
shear force, a^, to the outboard edge of all panels and the forward 
edge of the rib webs. From the equilibrium of the free body, forces 
a^, a^, and a^ may be found as functions of the unit force a^. 

1 b 



a^ = a, 1^ / 



a = a. W. / W 
C 11 c 



a, = a. 1 / W 

d id c 



( 20 ) 



The geometry of the panels and the load ratios are shown in Table III. 
3. Non-s 3 rrametrical shear panels. 

a) The cover skin panels have roughly the shape of a parallelogram, 
with the inclined sides skewed slightly. The maximum difference in 
inclination, cT- ^ = 3° 34*, occurs at the forward panels, elements 38 
and 15. These panels were treated in accordance with Reference 6 
which shows that the energy of the swept panel is given as 



££Cr - 



where = cot ^ f- cot S cot (0 cot ^6 
The forces on element 15 are shown in Figure 13. 



( 21 ) 

( 22 ) 




a. 






so that 

F 

n,n 

may be written for F 
/ 



Fig. 13 

FORCE OTAGRAM OF COVER PANEL 
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2 

of element 15. 
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20, 20 



So. S2S 



31 






b) The cover skin panels introduce additional forces due to interaction 
between torsion and bending. Twist of a swept wing produces bending and 
vice-versa. This interaction was introduced by considering the stress 
distribution in the rhomboidal cover skin panels, as shown in Figure 12. 
These panels were assumed to be parallelograms in order to simplify the 
analysis. Consider the plate shown, loaded with a uniform shear force, 
a , and axial load, f^. 

a ^ a 





fan A 




SWEEP COUPLING OF COVER PANEL 



The axial force, f^, was assumed concentrated at the flanges which were 
made up mostly of skin material, so that the cover panels carry only 
shear. However, it is seen in Figure 12 that the shear forces along 
the sides of the panel produce additional tensile stresses, f ^ > which 
do work on the displacements produced by the primary axial stress, f^. 
For this panel, the strain energy can be written. 



2£U = -f- + Ff^ ^ 

The third term in Equation 22 is due to the interaction between the 
axial and shearing forces and does not appear in an equation, ^wfitten 
for a rectangular panel. This term causes the stress distribution in 
this model to show a root stress concentration. Reference 3 states that 
’’Studies made on simple one cell swept box models with streamwise ribs, 
utilizing parallelogram shaped cover sheets, show that unless this cross- 
product term is introduced, there will be nothing in the analysis to 
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differentiate between a straight and swept wing, and there will be no 

3 

root stress concentration factor** . For the purpose of calculating the 
energies in the flanges adjacent to each panel, this axial load, f^ 

was split equally and added to the end of each flange. These coupling 
loads are shown as dashed vectors in Figure 6. 

C. Perform matrix operations. 

The matrix-force method, as described in References 1 and 3, 
expresses the energy of the entire structure as the successive product 
of three matrices. 

B£U = [a\ [f\ (a] ( 24 ) 

where a, a are the internal loads acting on the members. The 

In 

F matrix is symmetrical and represents the elasticity of the structure. 
The internal loads were found by writing the equilibrium equations for 
the determinate structure and for each redundant member. The number 
of equations must equal the total number of unknown internal loads less 
the number of redundants. The equations are written so that all the 
independent external loads and redundants are collected on one side and 
non-redundant loads on the other side of the equality; 
such as 

y] = [-S] i a,} (25) 

and 



where 



are matrices of the coefficients of the loads 



Solving for the non-redundant loads, a^^ , as functions of the applied 



loads, P , and the redundants, a^ , we get: 



(26) 
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Then, by adding the trivial relation between the redundants 



in 



( 27 ) 



(28) 



(29) 



and expanding ^jj to include Equation 4 results i 

= [<5] ('= i 

which may be transposed, yielding Equation 6 

{<2j - L^J “ I a,J [e]^ 

Now, substi'tuting into Equation 1, 

2£U = [P la^JsVs j 
= [«] 

is a symmetric matrix. Its members are the coefficients of the 
strain energy terms when the energy is written as a function of applied 
and redundant loads. Partitioning ^h| into four sub-matrices. 



(30) 




P p i 



H I 1-4 

rtp j rtr 



yP ^ appl/ed loads - 33 
r - tedunJani loads -/6 



Rewriting Equation 29, 

2.EU = 'h 2Pa^ •ha'l (31) 

Minimizing the strain energy with respect to the redundant, a^, yields 



Rearranging 





= n 2 


Z.OL H ^ 






/z r zz 




= — 




" K,] w 






£n\ 





(32) 



(33) 
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The matrix of all the unknown internal loads is given by 



[a] 



TOTAL 












( 34 ) 



Hor., - [S]W 



S is the unit load distribution matrix for the indeterminate structure. 

Each S. term is the value of the i*th member load when the structure 
im 

is subjected to a unit value of the m*th applied load. 

Using the notation of Reference 5, the strain energy can be 
written: 



c/ = i H [c] H 



(35) 



where C is the matrix of flexibility influence coefficients. Sub- 
stituting from Equation 11 into Equation 23, 



U = — 

2E 






therefore 



= [i- |£A/] 

- £ HMM 



_r 

^■/v' 



(36) 



(37) 



and substituting EN = — [H 22 j [^2l] Equation 10 

Note that Equation 34 may be expanded to read 



a = s, - P. + s. , P- . . . s, P 
1 11 1 12 2 1mm 



(38) 



(39) 



^2 ®21 ^1 ®22 ^2 • " ®2m ^m 
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Therefore, each internal load, a^ , may be found by substituting a 
particular set of applied loads, ^pj . For any other loading condition, 
only the matrix of applied loads needs to be changed to obtain the new 
internal loads for that loading (flight) condition. This enables a 
rapid evaluation of the maximum load on any element over a range of 
applied loads. 

A computer program, originally set up by Reference 2 for the CDC 1604 
computer was rewritten in Fortran IV to fit the IBM 360/67 computer and 
appears in Appendix A. This program performs the complete matrix oper- 
ation outlined in Section C with input A and B from the flexibility matrix 

»»» 

and the applied loads, P. Matrices A, .B, F, and P appear in Tables IV, 

V, VI, and VII. 
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EXPERIMENTAL RESULTS 



In order to ascertain whether or not the results obtained from the 
strain energy method were accurate for this type of thick skinned wing 
construction, it was necessary to have experimental data with which to 
compare the results of the idealized model. The wing is mounted in a 
jig which closely approximates the actual method of joining the wing to 
the aircraft fuselage, as shown in Figure 15. Since the wing is a 
continuous structure, the root fixity coefficient is not changed as when 
a single span wing is mounted on a completely rigid bulkhead. The 
loading condition investigated was a torque of 336,000 in. lbs. This 
load was considered to be well within the elastic range, but still of 
sufficient magnitude to give adequate strain levels throughout the 
structure for repeatable strain measurements. Some strain gages were 
installed on the wing when received from Ling-Temco-Vought Aerospace 
Corp., and others were installed at the Naval Postgraduate School. 

Figure 16 shows the wing with one cover skin removed for installation 
of strain gages. Wing section = 34.5 was particularly well instru- 
mented, with strain gage rosettes placed on both sides of top and bottom 
skins to check for differential bending effects. Strain gage data were 
reduced to principal stresses and directions through computer program 
**Recrose”, written in Fortran IV for the IBM 360/67 computer, as shown 
in Appendix B. The torque loads were applied with four hydraulic cylin- 
ders through a heavy fitting attached to the wing fold lugs. The applied 
torque was designed to be perpendicular to the elastic axis. The load 
was applied at both wing fold axes; and to ensure a S3nnmetrical load, 
the four hydraulic cylinders were supplied through a common manifold. 
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rtjy with Hover 3>in Ren^ved 



Load magnitude was determined by a manifold pressure gage in the 
hydraulic pressure system and checked by four Dillon dynamometers 
attached to each cylinder. The strain readings were taken with a 
Budd/Datran digital strain indicator (Model TC22). All gages were 
wired to a common junction panel as shown in Figure 17. All gages 
were balanced at a zero reading before applying the load. Eighteen 
gages were read at each loading in about one minute. This minimized 
the chance for drift in gage reading. All gages were rechecked for 
their zero reading after the load was released and sufficient time 
had elapsed for hysteresis effects to be minimized. All testing was 
done at night to eliminate the thermal effects of sunlight shining on 
the wing. Reference 2 contains all strain gage coordinates, gage fac- 
tors and resistances. It also shows that the strains are in the linear 
region by comparing the gage readings at a loading of 8000 psi (336,000 
in. lbs. of torque) with those taken at 7000 psi (294,000 in. lbs. of 
torque) and extrapolated to 8000 psi. 

Strain readings for the rosettes on the cover panels appear in 
Table VIII. The readings were averaged for the back-to-back gages 
and entered in program **Recrose" to compute the stresses. The normal 
and shearing stresses and the principal axes for the skin panels are 
also shown in Table VIII. 
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DISCUSSION OF RESULTS 



There are two areas of comparison in this study. The first is the 
comparison of the theoretical results with the experimental results, and 
the second is the comparison of the theoretical analysis of the wing 
with and without including the center section. 

In the first area, the stresses as found in Table IX show fairly 
close agreement. In the panels where the best stress agreement occurs, 
such as panel 38, the strain readings on both sides of the skin are 
nearly equal. The strain readings of back-to-back gages were averaged 
and used to find the experimental stresses at each station. Where the 
corresponding interior and exterior gages showed considerable variance 
in their strain readings, the experimental and theoretical stresses also 
had greater variance. The difference in strain readings is due to 
differential bending in the cover skin and possible torsional resistance 
of the pivot rib at the rear fitting, both of which were neglected in the 
theoretical analysis. The differential bending was neglected because 
the model used for the flexibility coefficient of the cover skin was a 
flat panel which carries only shear. The torsional resistance of the 
pivot rib was neglected because the effect of an unknown amount of slack 
movement about the wing hinge could not be incorporated into the 
theoretical analysis. The strain gages themselves are believed to be 
accurate to within 5 micro-inches of the true strain values. Strain 
readings found by Holgren and Comfort in 1963 were very close to those 
found in this study. The substitution of aluminum pads between the 
jig and the loading floor for the original pl 3 rwood pads, as recommended 
in Reference 2, does not seem to have made much difference in the strain 
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values. The original gages were carefully installed and waterproofed 
which seems to account for their longevity. Most of the gages have 



been in place for over four years. The normal expected life-span of 
this type gage in industrial use is one or two years. The first factor 
mentioned, that of differential bending and torsional resistance of the 
pivot rib, seems most likely to be the cause of the difference between 
theoretical and experimental results. The differences do not seem to 
be too great since some of the references mention that errors of 100% 
are sometimes found in the case of thin skinned structures. More con- 
clusive comparisons could be made if the torque loading were increased 
and larger strain readings were available. However, it was decided that 
a significant increase in the torque loading beyond the 336,000 in. lbs. 
used in this study might permanently warp the wing and render it useless 
for further study. 

The second area of comparison is between the theoretical results 

found in this study, and those of Reference 2. Including the center 

section, from = 0 to = 25.095 resulted in changing 25 flexibility 

influence coefficients and did not change the equilibrium equations. In 

effect, it moved the wing root from the pivot rib at Y = 25.095 to the 

w 

center rib at Y^ = 0. Most of the coefficients changed by 3 or 4%, and 
none of them changed by more than 13%; yet the matrix of resultant loads. 



designated as a^ (Reference 2). The greatest change was in the load 
on the rear intermediate beam, R.I.B. a^ (Reference 2) = 1010.4; and 
a^2 ” 10,263.5. This load is applied to element 24 and will influence 
the stresses on panels 23 and 25. The flexibility coefficient associated 




shows some considerable changes. In comparing ^ATj of Reference 
found in this study, the former loads will be 



2 with 
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with this force has a much smaller change, ^^ 3 ? ^3 (Ref. 2 ) = 13.3832 
and = 14.3406. The purpose of mentioning these changes is to 

show that the cumulative effect of minor changes in 25 out of 129 non- 
zero flexibility coefficients may be great. This means that small errors 
in deciding upon the true area and length of the elements of the ideal- 
ized model may be significant. When the actual structure to be idealized 
has irregular cut-outs and joints, and reduction in skin, web, and' 
flange areas occurs in steps rather than linearly, considerable effort 
will be involved in just constructing an idealized model. By examining 
Table X, it is seen that inclusion of the center section relieved the 
stresses at the leading edge of the wing which was the area of greatest 
discrepancy in Reference 2. Both normal and shear stresses are much 
closer to the experimental values. The stresses in the rear panels are 
also improved, but it is felt that the pins at point B and C in Figure 2 
fit snugly enough to offer some torsional resistance of the pivot rib at 
the rear fitting, which was not accounted for in the theoretical analysis. 
Tables XI and XII list the resultant loads, a^ through a^^ for the wing with 
and without center section included. The agreement between theoretical 
and experimental stresses throughout the wing seems to bear out the val- 
idity of the matrix-force method, and the approach as described in this 
study in the prediction of stresses in a wing prior to its construction. 

For greater accuracy, it would be necessary to use more elements which 
would require the calculation of additional equilibrium equations and hence 
expansion of the matrices. 
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CONCLUSIONS 



The matrix-force method of analysis using the strain energy theory, 
as proposed in Reference 1, is acceptable for a thick-skinned wing. 

It can be expected to predict the stresses in a complex, highly redun- 
dant structure, loaded in torsion, with a reasonable degree of accuracy. 
For increased accuracy, the long slender panels representing the cover 
skin should be broken up, and smaller elements used which would reduce 
the size of the structural grid but increase the number of unknowns. 
Reference 7 mentions a method of accomplishing this by dividing each 
panel into a number of narrow longitudinal strips which could be treated 
as rectangles with constant shear flow in each one. The difference in 
shear forces of adjacent strips would be taken by assuming concentrated 
stringer areas in the spanwise direction. This method would still 
neglect Poisson's effect. From the results of this study, it appears 
that a better panel would be one whose length and width were of similar 
size. 

In summation, it was found that the accuracy of the matrix-force 
method was improved considerably (compared with the results found in 
Reference 2) by the inclusion of 14 additional members which resulted 
in changes to 25 flexibility coefficients. Generally speaking, it must 
be possible to attain a very high degree of accuracy by fundamentally 
the same method as has been used in this study, by dividing the structure 
into a much larger number of elements. This involves a larger number of 
redundancies and much more numerical work. This could be programmed to 
be accomplished by the computer. An experienced stress analyst, who 
knows the stress distribution of most of the elements of the structure 
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on which he is working, would be able to concentrate the redundancies 
to the unknown points, thus reducing the computational work considerably. 
The experimental analysis would be improved by the addition of more strain 
gages on the inside of the wind cover skins. 

Because of the number of man - hours involved in removing and 
replacing the cover skins, however, it would not be practical to install 
additional gages. 
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TABLE III 



THE EDGE SHEAR FORCES ON THE PANELS 




Piece 

No. 


Wi 

in. 


Wc 

in. 


lb 

in. 


Id 
in* . 


II 


a.c/a-^= 

W^/Wq, 


aa/ai= 


1 


5.64 


6.82 


14. 3408 


14.3287 


2.1027 


.8269 


2.1009 


2 


6.82 


7.44 


13.5077 


13. 5041 


1.8155 


.9166 


1.8150 


3 


7.44 


7.66 


13.0307 


13.0303 


1.7011 


.9712 


1.7010 


4 


7.66 


7.54 


12.5756 


12.5754 


1.6678 


1.0159 


1.6678 


5 


7.54 


7.06 


12.1370 


12.1296 


1.7191 


1.0679 


1.7180 


6 


7.06 


6.30 


11.8046 


11.7984 


1.8737 


1. 1206 


1.87 27 


7 


5.64 


7.24 


79.3021 


79.3021 


10.9533 


.7790 


10.9533 


8 


6.82 


8.94 


74.9206 


74.9206 


8.3803 


.7628 


8.3803 


9 


7.44 


9.66 


71.3846 


71.3846 


7.3897 


.7701 


7.3897 


10 


7.66 


9.70 


68.5890 


68.5890 


7.0710 


.7896 


7.0710 


11 


7.54 


9.20 


66.4567 


66.4567 


7.2235 


.8195 


7.2235 


12 


7.06 


8.32 


64.9038 


64.9038 


7.8009 


.8485 


7.8009 


13 


6.30 


7.20 


63.8339 


63.8339 


8.8658 


.87 50 


8.8658 


15 


14.3287 


20.8705 


79.3021 


74. 9206 


3.7997 


.6865 


3.5897 


17 


13.5041 


19.0991 


74.9206 


71.3846 


3.9227 


.7070 


3.7375 


19 


13.0303 


17.7879 


71.3846 


68.5890 


4.0130 


.7325 


3.8559 


21 


12.5754 


16.5114 


68.5890 


66.4567 


4.1540 


.7616 


4.0249 


23 


12.1296 


15.2393 


66.4567 


64.9038 


4.3608 


.7959 


4.2589 


25 


11.7984 


14. 1418 


64.9038 


63.8339 


4.5894 


.8342 


4.5138 


28 


11.8046 


14.1748 


64.9038 


63.8339 


4.5788 


.8327 


4.5033 


30 


12.1320 


15.2456 


66.4567 


64.9038 


4.3590 


.7957 


4.2571 


32 


12.5756 


16.5132 


68.5890 


66.4567 


4.1535 


.7615 


4.0244 


34 


13.0307 


17.7879 


71.3846 


68.5890 


4.0130 


.7325 


3.8559 


36 


13.5077 


19. 1022 


74.9206 


71.3846 


3.9220 


.7071 


3.7369 


38 


14.3408 


20.8878 


79.3021 


74.9206 


3.7965 


.6865 


3.5868 
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MATRIX A INPUT FORM 
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O' 




mi^rgr- 


O' 


O' 


r*^rg<\j 


O' 


O' 



I I I 



^ rn sO rvj (\j Nt >t O o O oo fO >o h** o o O ^ ^ O ^ 
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fO t —4 C\l (\J (NJ OsJ t —4 f\J C\i C\J C\J 

r^m^^ijn^f^QCoocoooo^NJ^t >ocooooo^<Nj>t‘invor-‘Coa'a'a" c>oo^c\j(Njm 
^ ^ ^ ^ ^ (\j (\j fsj (\j (\j ovj rvj (\j cvj (\j rsj <\j <\j rvj rvj m ro m ro rn 
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^OlM^<^J<^JO^^OfOC^^o*-^co^t(\l^Oln^tt\lln ^mrsjrvjrvjmonj^o^ r- 
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1 1 

3 3 

5 5 

7 7 

8 4 

8 8 

812 

816 

820 

824 

828 

832 

836 

840 

844 

848 

945 

1123 

1237 

1349 

1538 

1642 

1829 

1934 

2010 

2014 

2018 

2022 

2026 

2033 

2039 

2043 

2047 

2121 

2125 

2129 

2133 

2142 

2146 

22 1 

22 5 

22 9 

2213 

2217 

2221 

2225 

2229 

2233 

2242 

2246 

23 B 

2312 

2316 

2323 

2324 

2328 

2332 

2341 

2345 

2349 

27 9 

29 3 

29 7 

2911 

2915 

2919 

2923 

2927 

29 31 

2935 

2939 

29 4 3 

2947 

3018 

3216 



TABLE V 

MATRIX B INPUT FORM 



1 .0 


134 


0 . 0 


2 20 


1.0 


244 


- 1.0 


1 


- 1.0 


335 


1.0 


4 4 


- 1.0 


436 


1.0 


1 


- 1.0 


537 


1.0 


6 6 


- 1.0 


638 


1.0 


1 


- 1.0 


8 1 


55.2466 


8 2 


69.5754 


8 3 


83.0796 


1 


96.1099 


8 5 


108.6853 


8 6 


120.8150 


8 7 


132.6135 


1 


7.4179 


8 9 


5.9191 


810 


4.2884 


811 


3.0187 


1 


1.8124 


813 


0.5448 


014 


0.5448 


815 


1.8124 


1 


3.0187 


817 


4.2884 


818 


5.9191 


819 


7.4172 


1 


2.5220 


821 


3.2013 


822 


3.5384 


823 


3.5707 


1 


3.3099 


825 


2.8758 


8 26 


2.3642 


827 


2.5220 


1 


3.2013 


829 


3.5384 


8 30 


3.5707 


831 


3.3099 


1 


2.8758 


833 


2.3642 


834 


- 15.9213 


835 


- 30.7843 


1 


- 45.6107 


837 


- 60.8687 


838 


- 75.9532 


839 


- 3.2013 


1 


- 3.5384 


841 


- 3.5707 


842 


- 3.3099 


843 


- 2.8758 


1 


- 2.5220 


845 


- 3.2013 


8 46 


- 3.5384 


847 


- 3.5707 


1 


- 3.3099 


849 


- 2.8758 


921 


1.0 


934 


8.3804 


1 


- 1.0 


1022 


1.0 


1035 


7.3897 


1046 


- 1.0 


1 


1 .0 


1136 


7.0710 


1147 


-1 .0 


1224 


1.0 


1 


7.2236 


1248 


- 1.0 


1325 


1.0 


1338 


7.8009 


1 


- 1.0 


14 7 


0.0 


14 8 


7. ''563 


1532 


- 1.0 


1 


- 7.8009 


1543 


1.0 


1631 


- 1.0 


1637 


- 7.2236 


1 


1 .0 


1730 


- 1.0 


1736 


- 7.0710 


1741 


1.0 


1 


- 1.0 


1835 


- 7.3897 


1840 


1.0 


1928 


- 1.0 


1 


- 8.3804 


1939 


1.0 


20 8 


- 0.99948 


20 9 


- 0.99980 


1 


- 0.99999 


2011 


- 0.99996 


2012 


- 0.99974 


201 3 


- 0.99915 


1 


0.99915 


2015 


0.99974 


2016 


0.99996 


201 7 


0. 99999 


1 


0.99980 


2019 


0.99948 


2020 


- 0.69666 


2021 


- 0.65009 


1 


- 0.60391 


2023 


- 0.55916 


2024 


- 0.51788 


2025 


- 0.48235 


1 


- 0.45373 


2027 


0.69666 


2028 


0.65009 


2029 


0.60391 


1 


0.55916 


2031 


0.51788 


2032 


0.48235 


2033 


0.45373 


1 


- 0.65009 


2040 


- 0.60391 


2041 


- 0.55916 


2042 


- 0.51788 


1 


- 0.48235 


2044 


0.69666 


2045 


0.65009 


2046 


0.60391 


1 


0.55916 


2048 


0.51788 


2049 


0.48235 


2120 


- 0.71732 


1 


- 0.75927 


2122 


- 0.79688 


2123 


- 0.82936 


2124 


- 0.85597 


1 


- 0.87645 


2126 


- 0.89114 


2127 


0.71732 


2128 


0.75927 


1 


0.79683 


2130 


0.82936 


2131 


0.85597 


2132 


0.87645 


1 


0.89114 


2139 


- 0.75927 


2140 


- 0.79688 


2141 


- 0.82936 


1 


- 0.85597 


2143 


- 0.87645 


2144 


0.71732 


2145 


0.75927 


1 


0.79688 


2147 


0.82936 


2148 


0.85597 


2149 


0.87645 


1 


56.885 


22 2 


56.885 


22 3 


56.885 


22 4 


56.885 


1 


56.885 


22 6 


56.885 


22 7 


56.885 


22 8 


1.83113 


1 


1.12519 


2210 


0.27134 


2211 


- 0.48011 


2212 


- 1.30551 


1 


- 2.34025 


2214 


- 2.34025 


2215 


- 1.30551 


2216 


- 0.48011 


1 


0.27134 


2218 


1.12519 


2219 


1.831 13 


2220 


2. 59681 


1 


3.39403 


2222 


3.84895 


2223 


4.02233 


2224 


3 . 9375 J 


1 


3.64622 


2226 


3.20813 


2227 


2.59681 


2228 


3.39403 


1 


3.84895 


2230 


4.02233 


2231 


3.93750 


2232 


3.64622 


1 


3.20813 


2239 


- 3.39403 


2240 


- 3.84895 


2241 


- 4.02233 


1 


- 3.93750 


2243 


- 3.64622 


2244 


- 2.59681 


2245 


- 3.39403 


1 


- 3.84895 


2247 


- 4.92233 


2248 


- 3.93750 


2249 


- 3.64622 


1 


56.8554 


23 9 


56.8736 


2310 


56.8844 


2311 


56.8827 


1 


56.8702 


2313 


56.8367 


2314 


- 56.8367 


2315 


- 56.0702 


1 


- 56.8827 


2317 


- 56.8844 


2318 


- 56.8736 


2319 


- 56.8554 


1 


0.0 


2321 


- 15.8462 


2322 


- 31.8510 


2323 


- 47.9019 


1 


- 63.5718 


2325 


- 78.4498 


2326 


- 92.3667 


2327 


0.0 


1 


15.8462 


2329 


31.8510 


2330 


47.9019 


2331 


63.5718 


1 


78.4498 


2333 


92.3667 


2339 


- 15.8462 


2340 


- 31.8510 


1 


- 47.9019 


2342 


- 63.5718 


2343 


- 78.4498 


2344 


0 . 0 


1 


15.3462 


2 346 


31.8510 


2347 


47.9019 


2348 


63.5718 


1 


78.4498 


2412 


6.8182 


2511 


7.4397 


2610 


7.6599 


1 


7.5385 


28 8 


1.0 


29 1 


1.0 


29 2 


1.0 


1 


1 .0 


29 4 


1 . :> 


29 5 


1.0 


29 6 


1.0 


1 


1 . D 


29 3 


0.03219 


29 9 


0.01978 


2910 


0.00477 


1 


- 0.00844 


2912 


- 0.02295 


2913 


- 0.04114 


2914 


- 0.04114 


1 


- 0.02295 


2916 


- 0.00844 


2917 


0.00477 


2910 


0.01978 


1 


0.03219 


2920 


0.01009 


2921 


0.01415 


2922 


0.01555 


1 


0.01487 


2924 


0 . )1249 


2925 


0.00971 


2926 


0.00705 


1 


0.01009 


2928 


0.01415 


2929 


0.01555 


2930 


0.01487 


1 


0 .01249 


2932 


0.00971 


2933 


0.00705 


2934 


- 0.76286 


1 


- 0 , 77''19 


2936 


- 0.78969 


2937 


- 0.81957 


2938 


- 0.84856 


1 


- 0.91415 


2940 


- 0 . J 1555 


2941 


- 0.01487 


2942 


- 0.01249 


1 


- 0.00971 


2944 


- 0.01009 


2945 


- 0.01415 


2946 


- 0.01555 


1 


- 0.01487 


2948 


- 0.01249 


2949 


- 0.00971 


30 9 


1.0 


1 


- 1.0 


3119 


1.0 


3117 


- 1.0 


3211 


1.0 


1 


- 1.0 


331 ? 


1 .0 


3315 


- 1.0 
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MATRIX F INPUT FORM 
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TABLE’ VII 



MATRIX ' ?• ' ■ 
Input Fbrm / . 



1 


1500.00 


2 


252,80 


3 


714.00 


4 


528.00 


5 


109.00 


6 


- 210.50 


7 


-27 82.00 


8 


- 4523.56 


9 


- 5212.17 


10 


- 5568.74 


11 


- 5395.86 


12 


- 4593.58 


13 


- 4654.94 


14 


- 4654.94 


15 


- 4593.58 


16 


- 5395.86 


17 


- 5568.74 


18 


- 5212.17 


19 


- 4523.56 


20 


2390.73 


21 


4825.58 


22 


5395.02 


23 


6125.17 


24 


6229.25 


25 


5802.52 


26 


2738.27 


27 


2390.73 


28 


4825.58 


29 


5395.02 


30 


6125.17 


31 


6229.25 


32 


5802.52 


33 


2738.27 



1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

2 
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TABLE VIII 










EXPERIMENTAL RESULTS AT 


v„ = 34. 


5 




PANEL 

NO. 


GAGE NO. 


STRAIN 

READINGS 


AVERAGE 


Cr max 




Inte- 

rior 


Exte- 

rior 


Inte- 

rior 


Exte- 

rior 


STRAINS 


O' min 


T^max 


15 




486 




61 


61 


1094.02 








485 




118 


118 




1008.22 






484 




-50 


-50 


-922.41 




21 


27 9 


459 


17 


25 


21 


851.68 






280 


458 


95 


100 


97.5 




726.87 




281 


457 


04 


-14 


-5 


-602.06 




25 


367 


468 


16 


-43 


-13.5 


678. 15 






368 


467 


-47 


-99 


-7 3 




627.44 




369 


466 


29 


11 


20 


-576.74 




28 


28 


474 


-19 


34 


7.5 


960.47 






29 


47 3 


137 


86 


111.5 




827.86 




30 


47 2 


27 


-8 


9.5 


-695.25 




30 


89 


477 


9 


9 


9 


457.66 






90 


47 6 


55 


63 


59 




512.27 




91 


475 


-12 


-20 


-16 


-566.87 




34 


19 


480 


-10 


-21 


-15.5 


599.43 






20 


47 9 


76 


65 


70.5 




568.22 




21 


47 8 


29 


10 


19.5 


-537.02 




38 


13 


483 


-80 


-76 


-7 8 


495.37 






14 


482 


61 


62 


61.5 




803.49 




15 


481 


35 


42 


38.5 


-1111.62 
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TABLE IX 

Comparison of Results at Station = 34.5. Loading 

336,000 in. lbs. of torque. Principal axis is measured 
from load reference axis, the Center Intermediate Beam, 
which has 33° 58' sweep back. 




MAX. or 


1 MIN. 


MAXIMUM SHEAR 






PANEL 

NO. 


normal! STRESS 


STRESS 


PRINCIPAL AXIS 




Experi- 


Theore- 


Experi- 


Theore- 


Experi- 


Theore- 




mental 


tical 


mental 


tical 


mental 


tical 


15 


1094.02 


1605.47 


1008.22 


1333.88 


31° 52' 


30° 12' 


17 




1530.10 




1092.17 




28° 16' 


19 




1047.07 




895.67 




38° 32' 


21 


851.68 


678.85 


726.87 


845.91 


CM 

in 

0 

o 


40° 43' 


23 




-792.63 




819.44 




48° 05' 


25 


-602.06 


-970.55 


627.44 


1033.12 


38° 48' 


49° 20' 


REAR BEAM 


28 


960.47 


767.42 


827.86 


870.20 


44° 44' 


47° 36' 


30 


457.66 


626.57 


512.27 


556.54 


39° 20' 


45° 24' 


32 




813.00 




669.19 




40° 12' 


34 


599.43 


851.95 


568.22 


761.82 


37° 26' 


39° 55' 


36 




-7 67 .60 




901.15 




35° 49' 


38 


-1111.62 


-1052.83 


803.49 


1160.03 


27° 11' 


33° 56' 
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TABLE X 

COMPARISON OF THEORETICAL RESULTS AT STATION = 34.5 

(1) WITHOUT INCLUDING CENTER SECTION, AS IN REF. 2. 

(2) INCLUDING CENTER SECTION 




MAXIMUM 


minimum 


MAXIMUM 


PANEL 


NORMAL 


STRESS 


NORMAL -STRESS ' - 


• SHEAR” 


STRESS 


NO. 


(1) 


(2) 


(1) 


(2) 


(1) 


(2) 


15 


2234 


1605 


-1280 


-1062 


1757 


1334 


17 


2025 


1530 


-651 


-654 


1338 


1092 


19 


1420 


1047 


-630 


-744 


1025 


896 


21 


1009 


67 9 


-67 9 


-1013 


844 


846 


23 


827 


846 


-619 


• -793 


723 


819 


25 


994 


1095 


-418 


-970 


706 


1033 


28 


813 


97 3 


-37 3 


-7 67 


593 


870 


30 


714 


626 


-478 


-486 


596 


556 


32 


835 


813 


-505 


-525 


670 


669 


34 


1068 


852 


-47 8 


-67 2 


773 


762 


36 


1575 


1034 


-575 


-7 68 


107 5 


901 


38 


1982 


1267 


-1004 


-1053 


1493 


1160 



56 



1 

2 

3 

4 

5 

6 

7 

8 

9 

1C 

11 

12 

1 ? 

14 

15 

If 

1 7 

18 

19 

20 

21 

22 

23 

24 

25 



TABLE XI 



matrix of all internal loads 

{AT}- [S] {P> 

FROM IDEAIIZED MODEL WITH CENTER SECTION INCLUDED 



174.54054260 



1325.45092773 


26 


266.98681641 


462. C38330C8 


27 


1781. 35253906 


474. 22583DC8 


28 


3944.20239258 


675 .57104492 


29 


-8407. 73046875 


i:66.64038G86 


3C 


-7190.99218750 


-597.00024414 


31 


-5754.24218750 


-3451.01147461 


32 


-4264.46093750 


-3498. 75415039 


33 


-2847. 58496094 


-4256. 8125COOC 


34 


650.99536133 


-5356.39062500 


35 


716.78784180 


-5181 . c39C6250 


36 


392.5336914 1 


-4882. 4R43750G 


37 


-1 36. 54725647 


-7350.36718750 


38 


1679.07592773 


-7291 . 8125000C 


39 


6767. 17968750 


-5829.66 796R75 


40 


6785.34375000 


-5258.5 1953125 


41 


5224.67187500 


-4711 . 71Q93750 


42 


985.44262695 


-4C 13.64619141 


43 


10263. 53515625 


- 3407.25048828 


44 


3625.59106445 


397':. 232421H8 


45 


6134.81640625 


2286.73657227 


46 


5103.66015625 


2' 93. 25366211 


47 


8034. 79296875 


-1755. 8449737C . 


48 


4810.60546875 


-57;.25122!77C 


49 


6020.46484375 
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1 

2 

3 

A 

5 

6 

7 

8 

<; 

1C 

1 1 

12 

13 

1 A 

15 

16 

17 

18 

IS 

2C 

21 

22 

2 3 

2-'h 

25 



TABLE m 



MATRIX OF ALL INTERNAL LOADS 
(AT}-[S]{P} 

FROM IDEALIZED MODEL WITH CENTER SECTION NOT 



-92.32380e7t 26 

l‘^92. 3 1 51 6 555 27 

-60.230A3823 28 

-222.629623A1 2S 

-2 53.088348?*^ 30 

-94.439C1C62 31 

442.25463867 32 

-2286.37548828 33 

-3785.3C151367 34 

-4119.C9765625 35 

-4763 .984375CC 36 

-4985.19921875 37 

-4969.66796875 38 

-4492 . 2734375C 39 

-435^.52343750 4C 

-488^.99609375 41 

-5013.0 1953 12*= 42 

-4765. 91C15625 43 

-42 19. 2539C625 44 

-3867 . ->4985352 45 

5143 .6875GCCC 46 

20C. 63256836 47 

-379.16625977 48 

-1720.40966797 49 

-1354.33C56641 



INCLUDED 

-llCC. 95019531 
-939.57446289 
-859.61694336 
-3524.74560547 
-3361.27685547 
-3227.84252930 
-2977.405C293O 
-2714.C4394531 
767.79956C55 
904.96679687 
592.9165C391 
-17.30769248 
-720.509C33?f 
8591.921875CC 
8726. 14Ct25CC 
7870.3671875C 
4969.57031250 
1C 1C .43847656 
546C . 734375CC 
8681 .625CCCCC 
a 3 5 7 . C 3 1 2 5 C r 0 
798C.82C3125C 
45C3. 76171875 
769.41943259 
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//ENERGY JOB A SZ 5 , 0 1 FP , 5 , 5 , MESSERSCHMIOT.O.fl. ,MSGLEVEL 1 
// FXEC FORTCLG 
//FORT.SYSIN DO • 

C COEFFICIENT MATRICES A, B, AND UC ARE BASIC STRUCTURAL DATA 

ODIMENSION A('5rj,50).AIf50,50),B(50,50),CJ(50,50'),UCC50,50'), ^ 
1GC50,50'),H ISO, 50).GT<'50,50),H22('50,50),H12<’50,56) ,UCG fSO.SO), 
2H21('50,50) ,H22 K50,50),ENa6,33) ,AR(16) ,UN I TN fSO , SO ) , S 1 50 , SO) , 

3 AT ISO), HI 2N ISO, SO) ,C ISO , SO ) , L I ST C50 ) , P 13 5) , ENP 116 , 33 ) 

OEOU I VALENCE lA, C J, UCG, H22 , H2 1 , S ) , I B , G, H 1 2 ,C ) , (UC ,H ) , 

1 (AI,GT,H22I .UN ITN.H12N) 

C CLEAR ARRAYS TO RECEIVE INPUT 

DO 10 1=1,50 
00 10 J=l,50 
An.J) = 0.0 
RH ,J ) - 0.0 
UC(I,J')= 0.0 
10 AI (I,J ) = Q.O 

C READ NON REDUNDANT MATRIX A 

31 CALL INPUT(50,A) 

GO TO (31 ,32'), NEXT 

C READ LOAD AND REDUNDANT MATRIX B 

32 CALL INPUTIS J, R") 

GO TO (32, 33), NEXT 

C FORM ENERGY COEFFICIENT MATRIX UC 

33 CALL INPUT(50,UC) 

GO TO (33,35), NEXT 

35 DO 36 I = 1 , A9 
DO 36 J = 1,A9 

36 UC(J,I) = UC(T ,J) 

C INVERT THE A MATRIX 

CALL GAUSS3(33,1.F-U,A,AI,KER) 

GO T0(365 ,36A),KER 
36A PRINT 151 

151 F0RMAT(24H MATRIX A IS SINGULAR) 

365 PRINT 200 

37 PRINT 38 

380FORMAT(64H A MATRIX INVERSE BY GAUSS3, 9 COLUMNS, IE 19, 10 18... 
IPER PAGE// 

PRINT 100,((AI( I,J),J=1,9) ,1=1,33) 

PRINT 20) 

PRINT 100,((AI (I,J),J-10, 18),I=1,33) 

PRINT 200 ^ 

PRINT 100,((AI (I,J),J^19,27),I=1,33) 

PRINT 200 

PRINT 3C0, C(aI ( I, J ), J=28,33), 1=1,33 ) 

DC FORMAT(/9F10.6 ) 

150 FORMAT(/11F10.6) 

200 FORMATdHl) 

C POST MULTIPLY AI WITH B MATRIX 

CALL MATMPY(33,33,A9,AI ,B,CJ) 

C FORM G MATRIX BY ADDITION OF TRIVIAL RELATIONS TO CJ 

DO 40 1=34,49 

40 CJ Cl , I) = 1.0 
DO 41 1=1 ,49 
DO 41 J=l,49 

41 G(I, J) = CJ(I, J) 

PRINT 200 

42 PRINT 43 

43 FORMATC50H G MATRIX .FORMED FROM CJ WITH TRIVIAL RELATIONS//) 
PRINT 100,((G(I,J ),J = 1,9),I=1,49) 

PRINT 200 

PRINT 100,((G(I,J),J = D, 18), 1=1,49) 

PRINT 200 

PRINT 100,(CG(I,J),J= 19, 27), 1=1,49) 

PRINT 20 J 

PRINT 100,((rG(I,J),J = 28,36), 1=1,49) 

PRINT 200 

PR I NT 100, (CGCI ,J), J = 37,4 5) , 1=1,49 ) 

PRINT 200 

PRINT 4 00, C(G (I , J ) , J=46,4 9 ), 1 = 1,49) 

300 format ( / 6F13 .7 ) 

^90 format (/ 4F13.7) 

C TRANSPOSE G 

on 44 1=1,49 
DO 44 J=l,49 

44 GTCI , J) = G(J, I ) 

C FORM H MATRIX BY TRIPLE PRODUCT (GTICUOCG) 

CALL MATMPY (49 ,49,49, UC,G, UCG) 

CALL MATMPY (49, 49,49, GT, UCG, H) 

PRINT 200 
PRINT 45 
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45 FnRMATf44H MATRIX H FORMED BY TRIPLE PRODUCT, GTXUCXG//) 
PRINT 5QT, (■(• HH ,J ) , J-1 ,9), I - 1, 49) 

PRINT 20D 

PRINT 50D,rrHn,J) ,J=1 0,18), 1-1,49) 

PRINT 200 

PRINT SOO.CT Hfl , J) , J= 19, 27) , I- 1,49) 

PRINT 20u 

PR I NT 500 , n ,J ) ,J- 28,36), I =1,49 ) 

PRINT 200 

PRINT 500 , TTHf I, J ),J = 37, 4 5), 1 = 1,49 ) 

PRINT 200 

PRINT 600 , fC Hn ,J ),J=46,49), 1 = 1,49') 

5DI' FORMAT ( /PF13. 4 ) 

6jo format ( / 4F13.4 ) 

C READ LOAD MATR IX P 

DO 510 1 = 1,35 

510 PM) = 0.0 

57 RBA0(5,58) I, (P(I)) ,NEXT 

58 FORMAT ( 12 , F12.4, II ) 

CO TO ( 57,581 ) ,NEXT 

C PARTITION H MATR IX 

581 DO 54 1=34,49 

DO 54 J=34,49 

I I = 1-33 

JJ = J-33 

54 H226I I , JJ ) = HTI , J ) 

PRINT 200 

521 PRINT 522 , ('('H22n , J) , J = 1 ,8) , 1 = 1 , 16 ) 

522 FORMATf /8F12.4 ) 

PRINT 2UO 

PRINT 5 22 , ( TH2 2 n. O'),.) 16), I = 1,16) 

C INVERT H22 MATRIX 

CALL GAUSS 3 (16,1.E-10, H22, H22I,KER) 

GO TOf525 , 524) ,KER 

524 PRINT 15j 
PR I NT 20 ' 

525 PRINT 522,('fH22I('I,J),J-l,8),I = l,16) 

PRINT 5 22, (fH22I (1 ,J ), J-9, 16) , I-l , 16) 

DO 53 1-34,49 

no 53 J=l,33 

KK = 1-33 

53 H2HKK, J) H(I ,J ) 

C FORM REDUNDANT LOAD COEFFICIENT MATRIX, CN. 

DO 540 1=1,16 
DO 540 K=l,33 
SUM = C . iJ 
DO 541 J-I,16 

54] SUM » SUM -h H22UI ,J)*H2KJ,K^ 

54r. Rf^PCItK:) = SU^ 

00 56 1=1,16 

on 56 J=l,33 

56 CN(I,J') = -1. • * ENP(I,J) 

C FORW MATRIX UNITN BY PLACING EN BELOW UNIT MATRIX 

UO 99 I = 1,33 
on 99 J = 1,33 
99 UM I TNH , J ) - ■> .'* 

on 62 I ^ 1,33 

62 UNITN (I , n - 1.3 

DO 63 1=34,49 

OU 63 J=l,33 

KK = 1-33 

63 UNITNfl , J) = FNCKK, J) 

C FORM s MATRIX . S = TOCUNITN) 

CALL MATMPY (49,49,33, G, UNITN, S) 

PRINT 20.’ 

PRINT 803 

8(’C format (43H STRESS COEFFICIENT MATRIX. S= (G)(UN1TN) //) 
PRINT 15- , (( S ( I , J) , J=1 , 1 1 ) , I = 1,49) 

PRINT ?0j 

PRINT150, ((S(I,J),J=12,22),I =1,49) 

PRINT 20.' 

PR I NT 150, ((S(I ,J),J=2 3,3^), 1=1,4 9 ) 

C RFOUNDANT LOAD MATRIX, AR 

DO I = 1,16 
SUM = f . 

OU 81 J = 1 , 33 

81 SUM = SUM + IN( I , J) *P (J) 

8 AK(I ) = SU*-’ 

DO 59 1=1,49 

59 L I ST ( I) = ) * I 
PRINT 21 ■ 
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PRINT 60 

600F0RMAT I50H REDUNDANT LOAD MATRIX. AR - H22I H21 P 
1 16H AR( 1-16)= ^( )) 

PRINT 700,(LISTCJ).AR(J) ,J=1,16) 

700 FORMAT (■ /13X, I2,F26.8) 

; TOTAL UNKNOWN INTERNAL LOADS. AT- fS)fP) 

DO 6A0 1=1, A9 
SUM - 0.0 
DO 6A1 J=l,33 

6Al SUM = SUM + S(I,J)*PfJ) 

64^ AT( I) = SUM 
PRINT 200 

PRINT 900 ,, , 

9OO0FORMAT ( 45H MATRIX OF ALL INTERNAL LOADS. AT- ( S K P ) // 

1 16H AT m -AC ) ) 

PRINT 1000, CLISTn ),ATn) , 1=1,49) 

I'^'OC FORMAT ( /13X,I2,F20.8) 
no 52 1=1,33 

DO 52 J=34,49 

52 H12(I,J-33) = H(1,J) 

E - 10300000.0 

INFLUENCE COEFFICIENT MATRIX. C - ( HU + fH12) CEN)) / E 

DO 90 I - 1,33 

DO 90 K = 1,33 

SUM - O.D 

DO 95 J = 1,16 

95 SUM = SUM i- H12H, J)*ENfJ,K) 

90 H12Nn,K ) = SUM 
SCALE = inooo. 

DO 70 1-1,33 

DO 70 J-1,33 

70 CCUJ) = (CH(I,J) + H12Nn , J ))/E)*SCALE 
PRINT 200 
PRINT 73 

730FORMAT (^5 H MATRIX OF DEFLECTION INFLUENCE COEFFICIENTS / 
135H C =CCH11 + (H12)(N ))/ E)*(SCALE)) 

PRINT 150, ((C(I,J),J= 1, 11) , 1=1, 33 ) 

PRINT 200 

PRINT 150, ((C(I,J),J= 12, 22), 1=1,33) 

PRINT 200 , 

PRINT 150, ( (C(l ,J) ,J=23, 33) , 1=1,33) 

RETURN 

END , . 

SUBROUTINE GAUSS 3 f N , E P , A S , X S , KER ) 

IMPLICIT REAL*8('A-H) ,REAL*8fO-P ) 

REAL*4 AS( 53,50), XS(50,5'^),EP 
DIMENSION(A 50 , 50) , X f50, 50 ) 

DO 1 1=1, N 
DO 1 J = 1,N 
A( I, J )=AS n , J) 

1 X ( 1 , J ) = 0 . 0 
DO 2 K=1,N 

2 X<'K,K)=1.C 
U 00 34 L=1,N 

KP=0 

Z=0.0 

DO 12 K=L,N 

IF (Z-DABS( A(K, L))) 11,12,12 

11 Z=DABS A K,L 
K P=^K 

12 CONTINUE 
IF(L-KP) 13,20, 20 

13 DO 14 J=L,N 
Z=A(L, J ) , 

A(L,J) = A('KP,J) 

U AfKP,J)=Z 
DO 15 J-1,N 
Z=X(L, J) 

X(L, J)=X(KP, J) 

15 X(KP,J) = Z 

20 IF fDABS('AfL,L)) -EP) 5 ), 50,30 

30 IF(L-N)31,34,34 

31 LPl-Li-1 

DC 36 K=LP1,N 
IF( A(K,D) 32,36,32 

32 KAT10=ArK, D/ACL, L) 

DO 33 J=LP1,N 

33 A(K, J)= A(K, J)-RATIO*AU, J) 

DO 35 J=1,N 

35 XCK,J)=XfK, J)-RATI0*XCL,J) 

36 CONTINUE 



// 
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CONTINUE 
^fO no AT 1= 1 , N 
I I-Nfl-1 
on A3 J=1 ,‘l 
S=^.T 

IFC 1 I-ri)Al,A3, A3 
A1 IIPl-Il+1 

DO A2 K= I I PI ,N 

A 2 s=s+ ACi 1 , K)*x(K, n 

A3 X Cl 1 , j ) »cxn 1 1 j)-s) /A(i 1 , 1 n 

KER=1 

I'O 51 1 = 1 ,N 

t>0 51 J=1,N 
51 XS (t , J ) = XCI . J ) 
kL TUKN 
AO KER=2 
RETURN 
END 

SLIPRCUTINF N>ATPPY (L, K, N, A, R, C) 
riMENSIOM AC-jT.S'tl, R(5'J,5,.0. CC5'-,501 
DO 1'. 1=1, L 
nn 1 . K=i ,N 

CUM = 3. j 

00 ?3 J=1,M 

2 ' SUM = SUM t A(1,J ) * eCJ,K) 

r c c I ,K ■) = suK 
KFTUPrj 
irjD 

SUBROUT INF ric>UT CN.A) 

UIMENSIOrj A{N,N),KC6) ^ 

Oroul VALENCE (KCD,Jl),fKC2),J2),(KC3),J3),CKfA),Il), 

1( K( 5) , 12 ) , (KC5 ) , I 3 ) 

1 READ C5,2n I1,J1,A1, I?,J2,A2, 13,J3,A3. 1A,JA,AA, NEXT 

21 format ( a (212, F12.61, 1 1 ) 

I'O n 1 = 1,6 

IF (_K( I ) .r.T. N 1 00 TO KO 
1'' CONTINUE 

IF ( I l.NF . '• ) A(I 1 , jn = A1 
IF C I2.NL. .n A(I2,J2) = A2 
IF CM.NF. O A(I3,J3 ) = A3 
IF ( lA.NE . .' ) A (lA, JA ) = AA 

OP TO (1,2), NEXT 

2 kFTURN 

I ' ■ PRINT 1 ) 1 , K , 1 1 , J 1 , A 1 , 1 2 , J 2 , A 2 , I 3 , J 3 , A 3 , I A , J A , AA 
1,1 FORMAT (' K=-,F5/ IX, Ar’I2,F12.6 )) 

STOP 
r NO 
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//RECROSE JUB ASZ5,01FP , MESSERSCHM 1 0 T , 0. 6 . ,MSGLEVEL 1 
// EXEC FORTCLG 
//FORT.SYSIN DO * 

RECTANGULAR STRAIN RUSETTE DATA REDUCTION FOR THE DETERMINATION 
OF PRINCIPAL AXFS, PRINCIPAL STRESSESt MAXIMUM SHEARING STRESSES, 
AND STRESSES ALONG THE ROSETTE AXES 1 AND 3. 



C 



/* 



READ C5i 100^ N, E,U 
L-0 
1«0 

1 I^I+l 

READ 5,101 IGAGE,E1,E2,E3 
E1>E1*1. jE- 06 
E2-F2»l .OE-06 
E3=E3*1. OF-06 

CORRECTION FOR A GAGE FACTOR = 1.95 

E1-E1*1 .03 

E2 = F2*1. :-3 

E3=E3*1 .03 

G=1 . )-U 

H-1.0 U 

A=E1+ F3 

B=E1-E3 

C = C2.i*E2)-A 

R = SORTfCB**21+CC**2)) 

TAUMAX=(E/(2.''*H1) *R 
SIGMAX=(E/2. .0 * CCA/G) t^R/H)) 

SIGMIN =CE/2. 0) *({A/G)- CR/H)1 
PHIPRR=U. B*CATAM2CC,P.')) 

IF CB.GE.0.0) GO TO 2 
PHIPKR=PH IPRR+1 .5707963 

2 PHIPRU=PHIPRR/ .I745329E-'^1 
D=T AUMAX*CUSf 2 . T*PH1PRR1 

S I GM A 1 = res I CM A X+ S I GM I N 1 / 2 . 0 

S IGMA3 = rcS IGMAXtS IGM INI /2.0) -D 
TAUO=-TAUMAX*S lNf2.0*PHIPRR) 

IF n.GT.l) GO TO 3 
WRITE C6,102) 

3 WRITE C6,n3) IGAGE,E1,SIGMAX,E2,SIGMIN,F3,TAUMAX,SIGMA1,PHIPRR, 
9SIGMA3,PHIPRD, TAUO 

L=Li-ll 

IF CL.LT.66l GO TO 4 
WRITE 16,102) 

L= 0 

4 IF CI.LT.N) '0 TO 1 

1-n format CIR,E17.7,F10.5) 

IJl FORMAT CI8,3F12.5; 

102 FORMAT (IHl) , , 

113 FORMAT CT6,'R(iSFTTF NUMBF R ' , I 8 , // , T 1 6 , El= ,F14.7, IN/IN ,T45, 

3 'SIGMAX=' ,F14.7, ' PSI' ,/,T16, 'E2= ',F14.7, 'IN /IN' ,T45, ’SIGMIN=' , 

3f 14.7, ’ PS I ',/ , T16, 'E3= ' , F14 .7, ' I N/ I N ' , T45 , ' TAUMAX= ' ,F14.7, 

3 ’ PSI ',//, Til, 'SIGMA U’,F14.7, • P S I’ , T45 , ' PH I PR=',F14.7,' RADIAN 

3S',/, Til, 'SIGMA 3=',F14.7, ' PS I ' , T45 , ‘ t>H I PR*' ,F14.7, 'DEGREES',/, 

3T13, 'TAU ',F 14.7, 'PSI ' ,// 

PFTUR \l 
END 
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ized model were found by using the IBM 360/67 computer to solve for the redun- 
dant loads. These stresses show favorable comparison with the experimental 
stresses determined from strain gages mounted on the cover skin of the wing. 
This study shows that the matrix-force method is applicable to thick-skinned 
wings . 
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